Effect of Ohmic environment on optimally controlled flux-biased phase qubit 
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We analyze the effect of environment on the gate operation of flux-biased phase qubits. We em- 
ploy the master equation for a reduced density matrix of the qubit system coupled to an Ohmic 
environment, described by the Caldeira-Leggett model. Numerically solving this equation, we eval- 
uate the gate error as a function of energy splitting between qubit states, junction capacitance, and 
temperature. The analysis is presented for single-quadrature microwave (control) pulses as well as 
for two-quadrature pulses, which lower the gate error significantly for idealized systems in the ab- 
sence of environment. Our results indicate that two-quadrature pulses outperform single quadrature 
pulses even in the presence of environment. 
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Superconducting circuits containing Josephson junc- 
tions are promising candidates for scalable quantum in- 
formation processing [J-Q. However, small separations 
between successive quantum energy states in these cir- 
cuits 0, @ do not permit selective manipulation of the 
qubit in a two dimensional subspace and results in a dy- 
namical leakage of a quantum state to a broader Hilbert 
space of the circuit [6|. To reduce this leakage, Motzoi 
et al. Ref. pfl proposed a Derivative Removal by Adia- 
batic Gate (DRAG) method, which reduces the gate error 
to 10 _ J for an experimentally optimal gate time of 6 ns. 
This error is well below the required error threshold of 
10 -3 for fault tolerant quantum computation [8]. 

In addition to the dynamic leakage, any realistic model 
of a qubit must also address coupling of the qubit to en- 
vironment, which leads to further destruction of qubit 
states. Several efforts have already been made towards 
the study of accurate control of a qubit system 0, 0, [l(| ■ 
However, the effect of an environment on optimally con- 
trolled qubit has only been studied in a phenomological 
model [7] , which leads to the evolution of density matrix 
of the qubit in Lindblad form [111 ]. 

In this paper, we resort to a microscopic approach 
to the modeling of the environment. We employ the 
Caldeira-Leggett model of the system-environment cou- 
pling 



12t I bil l to describe time evolution of a flux-biased 



phase qubit, driven by the DRAG pulses pj. Numerically 
solving equation of motion for the qubit density matrix, 
we study the dependence of the gate error on tempera- 
ture and environmental coupling strength. 

Model. A flux-biased phase qubit consists of a 
Josephson junction (JJ) embedded in a superconducting 
loop [2] . Finite resistance of the J J results in dissipation 
processes in the qubit and can be accounted for by the 
Caldeira-Leggett model [l2|, EH . The full Hamiltonian of 
the qubit and the environment is 



H = H q + P(t) + H n 



V. 



(1) 



The Hamiltonian of the qubit H q is written in terms of 
operators Q and 8, the charge and phase difference of the 



JJ respectively: 
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where L (C) is the loop inductance (junction capaci- 
tance), <j) ex t is the external magnetic flux applied to the 
phase qubit, Iq is the critical current of the JJ, and 
<f>Q = h/2e is the flux quantum. The qubit is capacitively 
coupled to microwave current source, used to induce co- 
herent transitions between the qubit states [2| . This cou- 
pling introduces time-dependent part in the Hamiltonian: 
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(3) 



Here I(t) = I x (t) cos u>dt + I y (t) smujdt is microwave cur- 
rent with frequency w<j. 

The environment is introduced as a set of har- 
monic oscillators (reservoir) with the Hamiltonian Hr = 

ELi( m «/ 2 ) (p„/ ra Q + Wq^)- Tne coupling between 
the qubit system and the reservoir is bilinear in the JJ 
phase S and oscillator displacements x a : 



N 

V = ^2 J a x a q 

a=l 
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(4) 



where parameters j a determine the coupling strength be- 
tween the qubit and reservoir mode a. 

Our goal is to describe the time evolution of the qubit 
density matrix p(t). The qubit is initially prepared in a 
pure state, corresponding to the density matrix /5(0). As- 
suming that the environment is in a thermal equilibrium 
at temperature T, the master equation for p(t) has the 
following form [14| : 



dp 
~dt 



1 



H q (t),p(t) -Ct{p}, 



(5) 
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and the dissipative term is 
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where is a Heisenberg operator. In Eq. ([5]), Jji(i) 
accounts for dissipative part of the dynamics and m(t) 
represents the quantum noise of the environment jl5| : 

/>oo 

rn(t) = h J{oj) [1 + 2 N(w)] cos wtdu , (7) 
Jo 



772 (i) = ih J(uj) shout dw 
Jo 



(8) 



The spectral density J(oj) = Yla=i 7a/ (2™ q w q )(5(oj— w a ) 
for an ohmic environment is 



J(w) 
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where £ is a dimensionless parameter, and u) s is a cutoff 
frequency that exceeds all other frequency scales of the 
system. The Planck's function N(u) = l/[exp(7kj/T) — 
1] defines an average excitation number of environment 
modes with frequency w. 



In experimental setup (16J, the "potential" part of H q 
in Eq. ([2]) has one deep minimum and another very shal- 
low minimum that disappears at the critical flux <f> c . Ex- 
ternal flux (f> cxt is chosen in such a way that only a few 
levels are localized in the shallow well, but these levels 
are still separated from levels localized in the deep well by 
impenetrable barrier (l7j |. As a result, we truncate the 
qubit Hamiltonian, Eqs. ([2]) and ([3]), to three localized 
levels and obtain the following Hamiltonian in energy- 
representation: 



H q (t) 



2 



H r , 



(10) 



where flj — is the projector for the j th level, aj = 

— 1| is the raising operator, a= (I x — U y )e %Udt /2 is 
the amplitude of microwave drive, Xj — <fio(j\5\j — 1} /2ttH 
is the matrix element of the phase operator, ujj = (ej+i — 
£o)/h, £j is an energy eigenvalue of time-independent 
Hamiltonian H q and H nr contains non- resonant terms. 
In this three-level model, the lower two energy levels 
comprise qubit space while the third level accounts for 
a leakage level. 

Gate error and DRAG method. In order to quantify the 
error during gate operation we use gate fidelity averaged 
over two initial input states in a two dimensional Hilbert 
space, similar to one defined in Ref. 18]: 



F„ 



Tr 



UidealPjU} deal Pj(tg) 



(11) 



Here 1) 1 ideal represents an ideal evolution, pj(t) is an ac- 
tual density matrix of a qubit system with pj(0) 



,(0) 



and p"p represents two initial axial states in a Bloch 
sphere. The gate error E is defined as E — 1 — F g . 

A simple approach to minimize leakage of quantum 
information from qubit subspace is to use a single- 
quadrature Gaussian envelope pulse given by 

/*(*) = h{t) = ie-C- 1 '/ 2 ) 2 ^ - B , I y (t) = , 

(12) 

where t g is a gate time and a = t g /2. For a NOT gate 
operation, which we choose to focus on without any loss 
of generality, constant A is defined by J Q a I„(t) dt = tt 
and B is chosen so that the Gaussian pulse starts and 
finishes off at zero. However, such a pulse shape still 
suffers a large gate error of about 10~ 2 as shown below. 

The DRAG method reduces the gate error to order 
of 10~ 5 for a gate time of 6 ns by using two quadra- 
tures and time-dependent detuning d\(t) = luq — ujd — 
(A 2 — 4)/ 2 (£)/4A, where the anharmonicity parameter 
A = lj\ — 2 ujq, and A measures relative strength of — > 1 
and 1 — » 2 transitions, that is, A = A2/A1. We note 
that the laboratory frame is more suitable for the solu- 
tion of the reduced density matrix of the qubit coupled 
to environment. We preserve the form of the quadrature 
amplitudes as in Ref. 0] for a Hamiltonian in rotating 
frame 
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Then, we obtain the following equation for the microwave 
driving frequency for the Hamiltonian Eq. ((3]) in the lab- 
oratory frame 

toj d (t) + u d (t) = uj - di(t) , uj d (0) = uj . (13b) 

Although the DRAG correction is successful in reduc- 
ing the gate error below the required threshold, a practi- 
cal implementation may not be feasible due to stringent 
requirement to vary microwave frequency. For this rea- 
son, we also consider two-quadrature pulses with fixed 
driving frequency uj d — loq [19]. To obtain relation be- 
tween I x ,y(t) components, we perform an adiabatic trans- 
formation D = exp [—iI x a/2A (&¥ + Act|)] of the rotat- 
ing frame Hamiltonian as in Ref. [7] along with a dimen- 
sionless parameter a, and a\ — — — cr^]. We then 
require that the imaginary parts of — > 1 and 1 — > 2 
matrix elements of the transformed Hamiltonian vanish, 
while 1 — > 1 matrix element is reduced to the second 
order in I x (t)/A and obtain: 



A 



A 2 
4 



(14) 



For t he p hase qubit, A ~ y/2, which implies a = 0.5 (see 
Ref. [20J for transmon qubits for which a = 0.4 owing 
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FIG. 1. (Colors online) (a) Gate Error vs. Alpha for Gate 
Times t g uio — 250 (dashed blue) and t g uuo = 350 (solid black), 
(b) Probability vs. Time for Temperature T = O.lhuuo and 
£ = 2. The qubit is initially prepared in |1) state (solid black), 
which relaxes to |0) state (dashed blue) due to dissipation. 



to different value of A). Similar result is obtained from 
direct numerical simulation of the gate error for different 
values of a. As shown in Fig. QJa), a minimum value 
of the error occurs at around a = 0.5 for gate times 
t g ojo — 250 (dashed blue) and t g ujo = 350 (solid black). 

Results. Numerical parameters used below in our simu- 
lation are chosen accordingly to the actual experimental 
setup: C = lpF, J = 1.5 /i A, fi L = 2-Kl L/(j>n = 3.2, 
and <t>ext = 0.955 4> c , where <p c is a critical flux [161 ]. For 
these parameters, we find luq ~ 39 GHz, A ss 1.41 and 
A « -2.4 GHz. 

In Fig. [21 we plotted the gate error for the DRAG 
pulses with and without time-dependent detuning. We 
find that pulses with two quadratures and fixed driving 
frequency (thin dashed red) perform much better than 
single quadrature Gaussian pulses (thin solid black), but 
are not as effective as pulses with double quadratures 
and time-dependent driving frequency (thin dashed-dot 
blue). 

In order to study the effect of dissipation on the DRAG 
pulses, we solve the master equation (0 numerically. 
First, we consider the relaxation of the qubit from the 
first excited state to the ground state in the absence of 
microwave drive, which is shown in Fig. [TJb). F° r this 
simulation, we choose parameter £ = 2 so that the re- 
laxation time T\ w 700 ns corresponds to experimentally 
observed decay time for phase qubits 



21|. We note that 



the spontaneous relaxation rate of the first excited state 
can also be evaluated from the master equation ([5]) 
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For the above choice of dimensionless coupling param- 
eter £ we study the effect of dissipation on the DRAG 
corrections. In Fig. [21 we observe a non-monotonic be- 
havior of the gate error with gate time for pulses with 
the DRAG corrections. We find that for shorter gate 
times, the DRAG correction with time-dependent driv- 
ing frequency is less affected by dissipation (thick dashed- 
dot blue). However, for longer gate times, dissipation 
has a substantial effect on two-quadrature pulses. For 
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FIG. 2. (Colors online) Gate Error vs. Gate Time with 
(thick lines) and without (thin lines) dissipation for Gaussian 
(a = 0.5 t g ) (solid black), Gaussian (a = 0.5 t g ) with first or- 
der DRAG correction and time-dependent driving frequency 
(dashed-dot blue), and Gaussian (a = 0.5 tg) with u>d — ivo 
and a = 0.5 (dashed red), all in the laboratory frame. 



tgUjQ = 250 (t g f=a 6 ns), the gate error increases from 
10 -5 to higher order of 10 -3 when dissipation is turned 
on for the same DRAG pulses with dynamical detun- 
ing. This increase in the gate error is due to the relax- 
ation of the qubit from the excited state to the ground 
state, which becomes prominent for longer gate times. 
For comparison, we plotted the gate error for three dif- 
ferent types of pulses: single quadrature Gaussian pulse 
(thick solid black), the Gaussian pulse with first order 
DRAG correction and time-dependent driving frequency 
(thick dashed-dot blue) and the Gaussian pulse with the 
DRAG correction and resonant driving frequency (thick 
dashed red). One can conclude from these plots that the 
performance of two-quadrature pulses without detuning 
is comparable to the DRAG pulses with dynamical de- 
tuning when dissipation is turned on. 

Next we study the effect of temperature on the gate 
error. The plot (Fig. [3]) for the gate error normalized 
around error at zero temperature reveals a monotonic 
increase of the gate error with temperature. The mono- 
tonic increase in the gate error is due to increase in the 
relaxation rate with temperature. In Fig. [3l we plotted 
the gate error for two different gate times: tgus = 150 
(dashed-dot blue) and t g u! = 350 (dashed black). 

We now compare the results of numerical solution of 
the master equation (J5J), and the simple picture of the 
error due to coupling to the environment in terms of 
the "Fermi- Golden rule" transition rates. Considering 
an environment at zero temperature and assuming that 
the contribution to the error E from the environment is 
small, E <C 1 , we can evaluate the error as the probability 
of an excitation of a reservoir mode during the qubit oper- 
ation, which happens with rate Y: E(T — 0) = Tt g pn(t), 
where pn(i) = / * 9 pu(t)dt/t g is the time-average of 
probability of qubit being in the first excited state. At 
finite temperature, the processes with excitation of en- 
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FIG. 3. (Colors online) Normalized Gate Error vs. Temper- 
ature for Gaussian (a = 0.5 t g ) with the DRAG correction 
and time-dependent driving frequency from the numerical 
simulation for gate times t g iuo — 150 (dashed-dot blue) and 
tgLdo = 350 (dashed black). Analytical rate equation estima- 
tion of the normalized gate error (solid red). 



vironment happen with rate T(T) = F[l + N(u)]. In 
addition, the qubit can absorb an excitation from the en- 
vironment with rate TN(u). Combining these processes, 
we obtain the following estimate for the gate error due 
to coupling to the environment: 



E(T) 

Ttn 



= [{1 + N(lu )} + A 2 JV(wi - w„)] Pu(t) 



(16) 



+ N{Lo )p o{t) 



For average occupation of the ground and the first excited 
states being « 1/2, and for a weak anharmonicity of the 
qubit system |A| <C loq, the gate error reduces to 



E(T) 
25(0) 



1 + 4iV(w ) 



(17) 



The estimated normalized gate error (solid red) is plotted 
in Fig. [3] together with the gate error obtained from nu- 
merical simulation. The rate equation estimation of the 
error is fairly close to the error obtained from direct nu- 
merical simulation for a longer gate time (dashed black). 
However, for a shorter gate time (dashed-dot blue), the 
estimated error deviates from the exact numerical sim- 
ulation considerably suggesting that the rate equation 
description may not be valid for shorter gate times and 
higher temperatures. 

Discussion and Conclusions. In comparison to single- 
quadrature pulses, two-quadrature microwave (control) 
pulses lead to significant suppression of the gate error. 
Despite the presence of dissipative environment, two- 
quadrature pulses reduce the gate error close to the de- 
sired threshold for fault tolerant quantum computation. 
At the same time, we also determined that optimal two- 
quadrature pulses with fixed microwave frequency pro- 
vide similar level of the gate error as two-quadrature 
pulses with time- varying frequency. 



In addition, we observed a monotonic increase of the 
gate error with temperature, which is due to increase in 
the relaxation rate with temperature. We found that 
temperature dependence of the gate error for longer 
pulses can be captured by a simple error estimation based 
on the rate equations. Nonetheless, the simple estima- 
tion of the error for shorter pulses differs from the gate 
error obtained from direct numerical solution of the re- 
duced density matrix. Hence, we concluded that full den- 
sity matrix solution is necessary to calculate the error for 
shorter gate times. 
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